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ABSTRACT 



Context. The outer Galaxy, where the environmental conditions are different from the solar neighbourhood is a laboratory in which 
it is possible to investigate the dependence of star formation process on the environmental parameters. 

Aims. We investigate the X-ray properties of NGC 1 893 , a young cluster (~ 1 - 2 Myr) in the outer part of the Galaxy ( galactic radius 
> 1 1 kpc) where we expect differences in the disk evolution and in the mass distribution of the stars, to explore the X-ray emission of 
its members and compare it with that of young stars in star forming regions near to the Sun. 

Methods. We analyze 5 deep Chandra ACIS-I observations with a total exposure time of 450 ks. Source events of the 1021 X-ray 
sources have been extracted with the IDL-based routine ACIS-Extract. Using spectral fitting and quantile analysis of X-ray spectra, 
we derive X-ray luminosities and compare the respective properties of Class II and Class III members. We also evaluate the variability 
of sources using the Kolmogorov-Smirnov test and we identify flares in the lightcurves. 

Results. The X-ray luminosity of NGC 1893 X-ray members is in the range 10 29 5 - 10 315 erg s _1 . Diskless stars are brighter in 
X-rays than disk-bearing stars, given the same bolometric luminosity. We found that 34% of the 1021 lightcurves appear variable and 
that they show 0.16 flare per source, on average. Comparing our results with those relative to the Orion Nebula Cluster, we find that, 
accounting for observational biases, the X-ray properties of NGC 1893 and the Orion ones are very similar. 

Conclusions. The X-ray properties in NGC 1 893 are not affected by the environment and the stellar population in the outer Galaxy 
may have the same coronal properties of nearby star forming regions. The X-ray luminosity properties and the X-ray luminosity 
function appear to be universal and can therefore be used for distance estimations and for determining stellar properties as already 
suggested by Feigelson and collaborators. 

Key words. Stars: coronae, pre-main sequence, luminosity function, open clusters and associations: individual : NGC 1893, X-rays: 
stars 



^ ' 1. Introduction 



x 



The formation of stars, the evolution of circumstellar disks and 
eventually the formation of planets are key research topics in 
modern astrophysics. So far, several star forming regions of dif- 
ferent ages and in different environment conditions have been 
studied in order to explore the parameter space as well as pos- 
sible. Although several questions have been answered, some as- 
pects of the star formation mechanism and initial stellar evolu- 
tion are still obscure. 

A key question is whether and how star formation de- 
pends on the environmental conditions. Since the Jeans mass 
dep ends on gas temperature , chemical abundances, and density 
rcf JElmegreeill2fMl200l . some dependencies on the physical 
conditions of star-forming clouds are expected. Moreover, the 
observed variation in the initial mass fun ction (IMF) at low stel- 
lar masses across d ifferent environments dRevle & Robinl,l2001t 
iGould et aTl[l9 98) suggests that gravity is not the only important 



parameter in the cloud fragmentation and subsequent evolution. 
As a consequence, other processes, e.g. turbulence a nd/or mag- 
netic field strength in the cloud have been considered dShu et all 



1987; Padoan & Nordlund, 1999: iMac Low & Klesserl ExM). 



* Tables [TJ and [3] are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



In this context, a natural way to explore the differences due to 
the environment is to look at the outer Galaxy. Molecular clouds 
and very young stellar associations are presen t also a t large dis- 
tances from the Galactic Center dSnell et all [2002), implying 
that star formation is occurring in these regions. Conditions in 
the outer Milky Way are thought to be less favorable to star for- 
mation, when compared with those in the solar neighborhood: 
the average surface and volume densities of atomic and molec- 
ular hydrogen are much smaller than in the solar neighborhood 
or in the inner Galaxy ( Wouterl oot et all 11 990), the interstel- 
lar radiation field is weaker dMathis et all 1 19831) . prominent spi- 
ral arms are lacking and there are fewer supernovae to act as 
external trigge rs of star formation. Meta l content is, on aver- 
age, smaller dWilson & MatteucciL E992), decreasing radiative 
cooling and therefore increasing cloud temperatures and conse- 
quently pressure support. Moreover, the pr essure of the inter - 
cloud medium in the outer Galaxy is smaller (Elmegreen, 1989). 
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Given these relevant differences, the study of star forming 
regions in the outer Galaxy is likely to be helpful to ascertain the 
role of environmental conditions in the star formation process. 
First of all, the environmenta l condition s could affect the evolu- 
tion of protoplanetary disks: lYasui et al.l d2009l) measured, in the 
extreme outer Galaxy, a disk frequency significantly lower than 
in the inner part of the Milky Way. The eventual dependence of 
disk lifetimes on metallicity could affect the formation of plan- 
ets and may give also insight into the "planet-metallicity corre- 
lation", since stars know n to harbor giant planets appear to pref- 
erenti ally be metal rich (ISantos et all 120031; iFischer & Valentil 
I2005L and references therein). M oreover, a recent theoretical 
works (lErcolan o & Clarke, 2009) shows that the disk lifetime in 
regions of different metallicity can give insight into the general 
mechanism of disk dispersal, allowing to discriminate between 
the two currently dominant models of disk dispersal: photoevap- 
oration and planet formation. 

In principle, the metal abundance differences between the in- 
ner and the outer Galaxy could also affect the coronal properties. 
Indeed, the thickness of the outer convection zone in low mass 
stars is a function of metallicity (Piz zolato et all 1200 ll) . This, 
in turn, likely changes the dynamo process itself, i.e. the X-ray 
emission engine. Moreover, differences in the composition of 
X-ray emitting plasma may produce other differences: X-ray ra- 
diative losses, for example, are dominated by radiation of highly 
ionized atoms and therefore we may expect changes in the emis- 
sion measure, as well as, in individual lines. Finally, if there was 
a difference in the disk fraction, we would expect an indirect 



ray luminosities than Class III stars (Stelzer & Neuhauser, 2001; 


Flaccomio et all 2003 £ 


ii Stassun et all 120041 iPreibisch et all 


2005: Flaccomio et al- 


2006: Telleschi etail [2007b thus a dif- 



ferent fraction of disked stars may affect the global X-ray lumi- 
nosity function. The assumption o f a universal X-ray luminos- 
ity function has been proposed by [Feigelson & Getrrianl d2005l) 
and used to estimate the distance (iKuhn et all 1201 Oh and total 



populations of sev eral young clusters (e.g. Getman et al., 2006; 
iBroos et all 120071) . In this context, the study of a cluster lying 
in a very different environment could be helpful in demonstrat- 
ing the universality of X-ray properties of young star forming 
regions and justify the use of X-ray luminosity functions to de- 
termine the properties of young clusters. 

In order to investigate the X-ray properties of a cluster in the 
outer Galaxy and compare them to similar cluster in the inner 
regions of the Galaxy, we have identified, as a suitable target, 
the young (~ 1-2 M yr) cluster NGC 1893, w hose galactocentric 
distance is 1 1 kpc dPrisinzano et all 1201 ll) . NGC 1893 is in a 
spiral arm different than that of the Sun and it is located near 
the edge of the Galaxy, where the environmental condition are 
quite different from those of the solar neighborhood. Situated 
in the Aur OB2 association toward the Galactic anti-center, 
NGC 1893 is associated with the HII r egion IC 410 and th e 
two pennant nebulae, Sim 19 and Sim 130 dGaze & Shainl[l952h . 
known as "Tadpole". It contains a group of early-type stars with 
some mol ecular clouds b ut onl y moderate extinction. S everal 
studies (IVallenari et all 119991 iMarco & Negueruefal [ 2002: 
Shar ma et all 120071 ICaramazza et all 120081 iPrisinzano et all 
201 1) demonstrated that star formation is still ongoing and there- 



fore we can study the X-ray properties of the PMS stars, with 
particular focus on the differences between the Class II and the 
Class III stars X-ray emissio n, given that t he dis k evolution may 
depend on the environment C Yasui et all 12009). The morphol- 
ogy, the age distribution of the clu ster and the star formatio n 
history have been studied in detail bv lSanz-Forcada et al .1 (1201 lh . 



They find that the cluster has not relaxed yet to a spherical dis- 
tribution, and it still has different episodes of stellar formation. 
The analysis of the age and disc frequency of the objects in 
NGC 1893, in relation with the massive stars and the nebulae, 
revealed an ongoing stellar formation close to the dark molec- 
ular cloud and the two s maller nebulae Sim 129 and Sim 130 
(ISanz-Forcada et ail 1201 lh . While both massive and low mass 
stars seem to form in the vicinity of the denser molecular cloud, 
Sim 130 a nd specially Sim 129 harb our the formation of low 
mass stars dSanz-Fo rcada et aTl 1201 lh . A parameter that is still 
v ery uncertain for this clu ster is the metallicity. As we discussed 
in Prisinz ano et alj d201 lh . the literature studies about the metal- 
licity of NGC 1893 have am biguous results . In a paper focused 
on the Galactic metallicity, Roll eston et alj d2000t) found for 8 
NGC 1 893 members a slight i ndication of under solar abun- 
dances, iDaflon & Cunhal d2004) show a marginal indication of 
subsolar metallicity for 2 members of the cluster, but these stud- 
ies do not show a strong evidence of subsolarme tallicity, there- 
fore, as assumed also in IPrisinzano et alJ (1201 lh . we adopt for 
NGC 1893 a solar metallicity. 

This study is part of a large project aimed at investigat- 
ing the star formation processes in the outer Galaxy and is 
based on a large observational campaign on the young clus- 
ter NGC 1893. The first results, concerning the joint Chandra- 
Spitzer large program "The Initial Mass Function in the Outer 
Galaxy: the star formi ng region NGC 1 8 93" (P I. G. Micela), 
have been presented in ICaramazza et al.l d2008l) . In that earlier 
work we have analyzed the Spitzer IRAC maps of NGC 1893 
joined to the Chandra ACIS-I long exposure observation, giving 
the most complete census of the cluster to date. The IRAC ob- 
servations of NGC 1893 allowed us to identify the members of 
the cluster with infrared excesses (Class 0/1 and Class II stars), 
while the ACIS observation was used to select the Class III ob- 
jects that have already lost their disks and have no prominent 
infrared excess. Although that first catalog of members was in- 
complete, the presence of 359 members indicates that NGC 1893 
is quite rich, with intense star forming activity despite the "un- 
favorable" environmen tal conditions in the outer Milky Way. In 
IPrisinzano et al.l (120 1 lh . we continued the study of the properties 
of NGC 1893 by using deep optical and JHK data and compil- 
ing a catalog extending from X-rays to NIR data. In this sec- 
ond work, we have assessed the membership status of each star, 
finding 415 diskless candidate members plus 1061 young stellar 
objects with a circumstellar disk or Class II candidate members, 
125 of which are also H a emitters. Moreover, optical and NIR 
photometric properties have been used to evaluate the cluster pa- 
rameters. Using the diskless candidate members, the cluster dis- 
tance has been found to be 3.6±0.2 kpc and the mean interstellar 
reddening E(B-V)=0.6+0. 1 with evidence of differential redden- 
ing across the whole surveyed region. These previous studies 
show that NGC 1893 contains a conspicuous population of pre- 
main sequence stars with a disk fraction of about 70%, similar 
to that found in clusters of similar age in the solar neighborhood. 
This demonstrates that, despite the expected unfavorable condi- 
tions for star formation, very rich young clusters can form also 
in the outer regions of our Galaxy. 

In the present work, as part of the same project, we present 
the in-depth analysis of 4 50 ks long Chandra AC IS-I observa- 
tion, already introduced in ICaramazza et all d2008h . focusing on 
the X-ray propertie s of the 1021 detected s ourc es in light of the 
previo us results of Caramazz a et alJ d2008l) and IPrisinzano et aTl 
(1201 lh . To put our work in context, we will compare the present 
study with the 850 ks long Chandra ACIS-I observation of the 
Orion Nebula Cluster (ONC), termed Chandra Orion Ultradeep 
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Project (COUP) (iGetman et alll2005h . The COUP with 1616 de- 
tected sources provided one of the most comprehensive dataset 
so far acquired on the X-ray emission of PMSs. Several studies 
were conducted on the 1616 COUP detected sources, in order to 
characterize the X-ray emission and to understand the corona l 
processes in the PMS stars. For example, Preib isch et al.l (2005) 
have analyzed the dependence of the X-ray luminosity on the 
stellar properties, findin g for the O rion stars that Lx is a func- 
tion of the mass, Wolk et al.l (120051) have studied the variability 
of the solar mass stars, in order to have in dication of the emis- 
sion of the young Sun, (Favata et an, l2005l) have investigated the 
geometry of flare loops, resulting in the first observational indi- 
cations of magnetic structure connecting the star and the disk. 
These are just few of the several results of the COUP published 
in more than 20 papers. Since COUP gives the most comprehen- 
sive view of magnetic activity in young stars ever achieved in 
the nearest rich cluster of very young stars, it is used here as a 
touchstone for comparison with NGC 1893. 

In the present paper, we analyse the properties of the X- 
ray sources in NGC1893, focusing, in particular, on the X-ray 
luminosity of low mass stars. We compare the X-ray proper- 
ties of the NGC 1893 members with infrared excesses ( here- 
af ter Class II stars) , as de fined in Caramaz za et al ] (120081) and 
in iPrisinzano et al.l d201 lb with the members without infrared 
excesses (hereafter Class III stars), that are classified as mem- 
bers of the cluster by means of their X-ra y emission (see 
ICaramazza et al 1 12008. 1 Prisinzano et al lHoIl for details). The 
outline of the present paper is the following: in Sj2]we describe 
the data reduction, the photon extraction procedure, and present 
our X-ray catalog. In |j3]we describe the spectral analysis that 
we pursued by means of spectral fitting and quantile analysis of 
X-ray spectra. In |j4]we compare the X-ray properties of Class II 
and Class III of NGC 1893 and in Sj5]we present the variability 
properties of the members. In Sj6]we discuss our results and in S]7] 
we summarize our work. 

2. Data Reduction 

The X-ray observations of NGC 1893 combine four nearly con- 
secutive exposures of the cluster taken in 2006 November and a 
fifth exposure taken in 2007 January for a total exposure time of 
~ 440 ks. The combined X- ray image is shown in Fi g. Q] Source 
detection was performed by ICaramazza et al] (120081) and we re- 
fer to that work for details of the observations and the detection 
method. 

Starting from the ICaramazza et al ] (120081) catalog, we pro- 
ceeded to the photon extraction of t he source photons b y means 
of ACIS Extract.] (AE) v3.131 dBroos et all 120021) . an IDL 
based package of tools that can assist the observer in performing 
the many tasks involved in analyzing a large number of point 
sources observed with the ACIS instrument on Chandra. AE 
makes extensive use of TARA 0, a package of tools in the IDL 
language for visualizing and analyzing X -ray astronomical data 
in FITS format , CIAO dFruscione et all 12006) Fl the data anal- 
ysis system written for the needs of users of the Chandra X-ray 
Observatory , and FTOOLSQ, a general package of software to 
manipulate FITS files. 

The extraction of point sources performed with AE takes 
into account the Point Spread Function (PSF) for single sources, 

1 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 

2 http://www.astro.psu.edu/xray/docs/TARA/ 

3 http://cxc.harvard.edu/ciao/index.html 

4 http://heasarc.gsfc.nasa.gov/docs/software/ftools/ 



which strongly depends on the off-axis distance (ff). While the 
PSF is narrow and approximately circular in the inner part of 
the field of view (9<5'), it has a non-Gaussian shape at large 
off-axis, becoming broader and more asymmetric. 

Moreover AE calculates the shape of the model PSF at each 
source position (by means of the CIAO task MKPSF), consider- 
ing also that the fraction of the PSF where photons are extracted 
depends on crowding. In this way, AE finds a compromise be- 
tween having the largest PSF fraction (thus providing good pho- 
ton statistics for further spectral and timing analysis) and having 
the best signal-to-noise ratio, not to mention the importance of 
avoiding the contamination from nearby sources. 

After the calculation of the PSF shape, AE refines the initial 
source positions, in ou r case originally estimated by PWDetect 
(iDamiani et all Il997al) assuming a symmetric PSF, by correlat- 
ing the source images with the model of local PSFs. Following 
AE science hintfl this last procedure was only used for those 
sources lying off-axis by greater than 5 arcmin (321 sources), 
while for the rest of the sources we simply adopt mean photon 
positions. Coordinates listed in Table [TJ as well as their lcr un- 
certainties, are the result of this process. 

After re-computing positions, AE defines source extraction 
regions as polygonal contours of the model PSF containing a 
specified fraction of source events (fpsF)- Generally, we chose 
fpsF=90%, and computed the contours from the PSF for a mono- 
energetic source with E = 1 .49 keV. For a few sources in the 
denser parts of the field of view this fraction was reduced so 
as to avoid contamination with other nearby sources, in the most 
extreme cases down to fpsF ~40%. 

Although the ACIS-I instrumental background level is spa- 
tially quite uniform, the actual observed background varies sub- 
stantially across the NGC 1893 field due to the extended PSF 
wings of bright sources and to their readout trails. Background 
was therefore estimated locally for each source, adopting once 
again the automated procedure implemented in AE, which de- 
fines background extraction regions as circular annuli with in- 
ner radii 1 . 1 times the maximum distance between the source 
and the 99% PSF contour, and outer radii defined so that the re- 
gions contains more than 100 "background" events. In order to 
exclude contamination of the regions by nearby sources, back- 
ground events are defined from an image that excludes events 
within the inner annular radii of all the 1021 sources. 

Results of the photon extraction procedure are listed in 
columns 7-10 of Table Q] where we give the background- 
corrected extracted source counts for 0.5-8.0 keV, the associated 
error and the total background counts expected. In summary, our 
1021 X-ray sources span a wide range of photon flux, from ~3 to 
~7100 photons during the exposure time. Most sources are faint 
(e.g. 60% have fewer than 50 photons). 

In the last six columns of Table Q] we list some char- 
acteristics of the sources: the significance (signal-to-noise ra- 
tio), Kolmogorov-Smirnov probability that the lightcurves of the 
sources are constant, the median photon energy, the absorption- 
corrected X-ray luminosities and a flag indicating the method 
used to determine them. 



3. Spectral Analysis 

In order to determine the intrinsic luminosity of our sources, 
spectral analysis is necessary. Unfortunately, spectral fitting can 
be applied only to high counts sources while relatively faint 



5 http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide/node35 
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Note 1. Tahlefflis available in its entirety in the electronic edition. 
Col. (1): X-ray catalog sequence number, sorted by RA. 
Col. (2): IAU designation. 

Cols. (3) and (4): Right ascension and declination for epoch (J2000.0). 

Col. (5): Estimated standard deviation of the random component of the position error, ycr^ + cr*. The single-axis position errors, a~% and cry, are estimated from the single-axis standard deviations of the PSF inside the extraction region and the number of counts extracted 
Col. (6): Off-axis angle. 

Cols. (7) and (8): Net counts extracted in the total energy band (0.5-8 keV) in the extraction region; average of the upper and lower Icr errors on col. (7). 
Col. (9): Background counts expected in the source extraction region (total band). 
Col. (10): Net counts extracted in the hard energy band (2-8 keV). 

Col. (10): Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. A reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded region. 
Col. (1 1): Photometric significance computed as net counts divided by the upper error on net counts. 
Col. (12): Kolmogorov-Smirnov probability that the lightcurves of the sources are constant. 

Col. (13): Effective exposure time: approximate time the source would have to be observed on-axis (no telescope vignetting) on a nominal region of the detector (no dithering over insensitive regions of the detector) to obtain the reported number of counts. 
Col. (14): Background-corrected median photon energy (total band). 
Col. (15): Absorption-corrected X-ray luminosities. 

Col. (16): Origin of the X-ray luminosity, in col. (15): (a) from fit to a one temperature model spectrum; (b) from analysis of quantiles; (c) from a conversion factor derived from the analysis of the quantiles. 
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Fig. 1. True-color image of the young cluster NGC 1893. A kernel smoothing has been applied to highlight point sources. The 
energy bands of the RGB image are 0.5-1.1, 1.1-1.9, and 1.9-8.0 kEv for red, green and blue colors. Brightness is scaled to the 
logarithm of the photon number in the displayed pixel. The color model depicts zero flux as black. 

sources with poor statistics cannot be investigated using this 
method. 

To overcome this limit, we decided to perform spectral fit- 
ting C ^3. lb on 311 sources for which it is possible to bin the 
spectrum so that, in the background-subtracted (net) spectrum, 
each rebinned channel achieves a signal-to-noise ratio at least 3. 
In situations where the background level is high, this approach 
can produce higher quality binning compared to fixing a mini- 
mum number of counts in the source spectrum. For less intense 
sources we analyzed the spectrum of sources considering the 
quantiles of the distribution of energy for source counts. As we 
will describe in the follow ing for sources do wn to 30 counts, ap- 
plying a quantile method dHong et al.Ll2004l) . we have been able 
to derive the values of Nh and kT interpolating directly from 
a thermal grid of models with a good accuracy, while for very 
low statistics sources we have used the quantile analysis to de- 
rive a median conversion factor from the count rate to the X-ray 
luminosity. Energy quantiles have been also used to discern the 
extragalactic contaminants. 



3.1. Spectral Fitting 



Reduced source and background spectra in the 0.5-8.0 keV band 
have been produced with AE (see fj2|, along with individual "re- 
distribution matrices files" (RMF) and "ancillary response files" 
(ARF). 



We fit our spectra assuming emission by a thermal 
plasma, in collisio nal ionizati o n equ ilibrium, as modeled by 
the APEC code ( Smit h et all [2001). Elemental abundances 
are not easily constrained with low-statisti c s spec tra and were 
fixed at Z=0.3 Z (see Prisinzano et al. I Hn Ih, with solar 
abundance ratios taken from An ders & Grevessd (1989). The 
choice of sub-solar abundances is sugge sted by several X- 
ray studies of star f orming regions (e.g. Feigel son et al.L [2002; 
Magg io et all |2007). Absorption was accounted for using the 
WA BS model, parametrized by the hydrogen column density, 
N H dMorrison & McCammor3,ll983h . 

We fit source spectra with one-temperature (IT) plasma 
models using an automated procedure. In order to reduce the 
risk of finding a relative minimum in the x 1 spaces, our proce- 
dure chooses the best fit among several obtained starting from 
a grid of initial values of the model parameters: log(NH) = 21.0, 
22.0 cm" 2 and kT = 0.5, 1.0,2.0, 10.0 keV. 

Among the 311 sources for which we conducted the spec- 
tral fitting, for 286 sources it was possible to determine Nh and 
kT. Using the spectral fit parameters, we have then derived X- 
ray luminosities for each of the 286 sources modeled with one- 
temperatu re spectrum, consideri ng for NGC 1893 a distance of 
3.6 kpc (Prisinzano et all 1201 lh . This sources are tagged with 
an "a" in column (16) in Tableland their luminosity is listed in 
column (15). The 25 sources for which we were unable to fit the 
spectra with a one temperature model, have been analysed with 
the study of the energy quantiles of the spectrum (see following 
sections). 
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Median Energy Distribution 




Source Median Energy[keV] 

Fig. 2. Cumulative distribution function of the photon median 
energy of each source computed in the [0.5-8.0] keV range for 
all the sources. 



3.2. Quantile Analysis 

Spectral fitting analysis cannot be applied to faint sources, but a 
rough analysis of the whole sample can be achieved by the study 
of the median photon energy of each source. We plot in Fig.|2]the 
distribution of the source median energy computed in the [0.5- 
8.0] keV range. From this plot, we can retrieve some information 
about the population of our X-ray sources: the median value of 
the distribution is 1 .6 keV, the bulk of the distribution is between 
1 and 2 keV (62% of the median photon energy of the sources 
lies in this range), while 10% of the values are over 4 keV; this 
energetic tail of the distribution is probably due to extragalactic 
contaminants, even if we cannot exclude at this stage that stellar 
flares on NGC 1893 members play some role. While this kind of 
analysis is very helpful to understand the behavior of the entire 
sample, it does not help much for the characterization of single 
sources. 

The conventional spectral classification of weak sources is 
often performed in the literature by means of the ratios of source 
counts in different spectral bands, i.e. X-ray h ardness ratios 
(XHR) which can be consid ered "X-ray colors" dSchulz et all 
1989: lPrestwich et all 120031) . The flaw in this method is that the 
choice of the bands is strongly related to the spectral shape and 
that, depending on the number of counts in each band, the error 
bar related to the ratio may be very large. For these r easons, we 
decide to apply an alternative method, introduced by Hon g et all 
(12004 . that uses the energy value that divides photons into pre- 
determined fracti on, instead of th e ratio of counts in prefixed 
bands. Following Hong et al. ( 2004J), if E x % is the energy below 
which the net counts are x% of the total counts, quantile Q x is 
Qx = E f % Zf" i where Ei and E u „ are the lower and upper bound- 
ary of the full energy band, i.e. 0.5 and 8.0 keV in our case. The 
fractions we choose are the median (50% quantile), and the quar- 
tiles (25% and 75% quanti les). Since quantile s, for a given spec- 
trum, are not independent, iHong et aD (120041) selected the inde- 
pendent variables log{Qsol{\ - Q50)) and 3(225/275), m order 
to have information on the photon population in several portions 
of the spectrum. The derived quantile values will be compared 



Median Energy [keV] 
13 3 4 




-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

i°g(Q»,/(i-Q s o)) 



Fig. 3. Quantiles of X-ray spectra. The dots indicate the quan- 
tities derived from source energy quantiles. The two grids 
shows the same quantities derived from simulated spectra. 
The solid line grid refer to thermal spectra with Nh - 
[10 20 , 10 21 ,0.6 • 10 22 , 1.1 • 10 22 ,3.4 • 10 22 ,6 • 10 22 ] cot" 2 and 
kT = [0.3,0.5,1.1,2.3,4.0, 10.0] keV, while the dashed line grid 
to simulated power law spectra, typical of AGN, with the same 
Nh values and power law indices F = [0.0,0.4, 1.0, 1.6,2.0]. 



with those calculated from simulated spectra in order to derive 
the best spectral parameters. Figure [3] shows the scatter plot of 
quantiles measured for our sources, compared with grids calcu- 
lated from simulated spectra. The solid line grid refers to thermal 
spectra, while the dashed grid to power law sp ectra, with spec- 
tral index typical of AGNs (Bran dt et all 1200 II) . We note that, as 
expected, most of the sources lies in the thermal grid region. By 
interpolating the loci of our sources inside the thermal grid it is 
possible to attribute values of Nh and kT to each source and de- 
rive the intrinsic X-ray luminosity at the distance of NGC 1893. 
In figure|4]we compare, for the brightest sources, the X-ray lumi- 
nosity derived from the quantile analysis with that derived from 
spectral fitting. The luminosities derived spectral fitting are gen- 
erally consistent with those derived from quantile analysis and a 
linear regression between the two luminosities results in a coef- 
ficient 1 .04 ± 0.4 giving strong support to the effectiveness of the 
method. The sources for which the luminosity was derived from 
the quantile analysis are indicated with a "b" in column (16) of 
tabled] 

By analysing the sources outside the region covered by the 
thermal grid, we found three subclasses of objects: 

- Two temperatures spectra. There is a small fraction of fairly 
strong sources (between 30 and 50 photons) that lie over the 
left bound of the thermal grid. Due to the concavity of the 
grid, this is the locus where thermal spectra characterized by 
two temperatures lie. For these cases, in our one temperature 
approximation, we associated them to the lowest Nh curve, 
interpolating just the kT value, that depends basically from 
the median energy of the spectrum (see Fig. [3]). 

- Very poor statistics spectra. Most of the sources outside the 
thermal grid are sources with less than 30 photons, with 
very large errors on quantiles. For these objects, we decided 
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Fig. 4. Scatter plot of the X-ray luminsity derived from spectral 
fitting versus that derived from quantile analysis. The solid line 
indicates equal values of Lx derived with the two methods. 

to adopt a single median conversion factor from count rate 
to luminosity. The conversion factor was derived from the 
quantile analysis of sources with more than 30 photons and 
applied to sources with less than 30 photons. These sources 
are indicated with a "c" in column (16) of TableQ] 
- Power law spectra. A small fraction of the sources lie within 
the power law grid and will be discussed in £ 13.2.11 

3.2.1. Rejecting extragalactic contaminants 

As shown above, with the energy quantiles analysis it is possible 
to attribute a luminosity to each source in the sample. It is also 
useful to investigate the contaminant population in the sample. 
NGC 1893 lies toward the galactic anticenter, therefore we ex- 
pect to detect a large number of extragalactic sources. Moreover, 
the long exposure necessary to better investigate the faint pop- 
ulation of our distant cluster results in a high probability of de- 
tecting AGNs. In order to estimate the number of AGNs, we con- 
sidered the sensitivity for three different off axis regions in the 
Chandra ACIS-I field of view and calculated the minimum flux 
that we expect from a dete cted extragalactic source with a typical 
power law index F = 1.4 (Br andt et al.Ll2001|) in the direction o f 
NGC 1893, (N H = 6.56 ■ 10 21 cm' 2 . iDickev & LockmarAll990h . 
By comparing the obtained fluxes wi th the log N-Log S obtained 
from the Chandra Deep Field North dBauer et all 120041) . we ex- 
pect to find ~ 190 AGNs among our sources. 

From the diagram in Fig. [3] we find 128 sources that are 
outside the thermal grid and seem to be compatible with the 
power law grid. We have tagged them as candidate AGNs in 
our catalog, but we have not excluded them from further analy- 
sis. Caution is needed because a very hard spectrum, due e.g. to 
flares can be confused with a power law spectrum; indeed, from 
the analysis of the variability, we have found that 16 sources 
tagged as candidate AGN show big flares. The spectral analysis 
finds support from the spatial distribution of candidate AGNs. 
Looking at Fig. Q] the morphology of the cluster is evident and 
it is clear that the distribution of sources in the four instruments 
chips is not uniform. We compared the number of sources de- 



tected in each chip with the number of candidate AGNs. As we 
expected, the total number of sources is higher in chips II and 13, 
where the peak density of the cluster lies, while the spatial distri- 
bution of the candidate AGNs is different and not related to the 
spatial distribution of the cluster. Again, we note that the candi- 
date AGNs that show flares follow the same space distribution of 
the cluster, therefore we treat this population with caution but we 
do not exclude it from the sample. This analysis is summarized 
in Table |2 

3.3. Upper limit to the X-ray luminosity of undetected 
NGC 1893 members 

The purpose of the present work is to examine the X-ray proper- 
ties of the NGC 1893 members, investigate the properties of stars 
in different evolutionary stages and compare them with those of 
other star formation regions. For these reasons , we ta ke advan- 
tage of the work described in IPrisinzano et alj d201 lb in which 
the optical and the infrared properties of NGC 1893 members are 
described. In order to compare the luminosity of stars from dif- 
ferent classes without any bias due to different depth of infrared, 
X-ray and optical images, we derived u pper limits to the X-ray 
luminosity for all the cluster members in Prisinza no et alJ d201 lh 
falling in the Chandra ACIS-I field of view and undetected in X- 
rays. Uppe r limits to the photon co unt rates were calculated with 
PWDetect dPamiani et all Il997allbh . with a detection threshold 
significan ce set to 4.6cr , the same used for source detection 
dCaramazza et all 120081) . In order to convert count rate upper 
limits into X-ray luminosity for X-ray undetected members, we 
calculated a conversion factor, considering the median value of 
the ratio between X-ray luminosity and count rate for the de- 
tected X-ray sources as done for the detected source s. In table 
m we l ist the sequence number of the source in IPrisinzano et all 
(201 1), the coordinates and the upper limit to the X-ray luminos- 
ity of the 158 Class II members that we will analyse. 

4. X-ray properties of low mass members 

Focusing on the X-ray properties of the NGC 1893 low- 
mass members, we started from the member catalog of 
Prisi nzano et alJ d201 ll) . In particular, we selected all the Class 
II or Class III candidate members with mass smaller than 2 M Q 
and whose optical colors are compatible with the cluster locus. 
In order to have a complete sample of sources, we selected all 
the stars with mass greater than 0.35 M Q , the choice of this com- 
pleteness threshol d derives from analys is of the color-magnitude 
diagram (fig. 1 1 in lPrisinzano et all 1201 ll) and will be confirmed 
in the following. We also required that the candidate members 
are inside the Chandra ACIS-I field of view. In this way, we 
identify 591 low-mass members, 307 of which are Class II stars, 
while 284 are classified as Class III stars. Among the 591 mem- 
bers, 158 of the Class II stars are not detected in the X-ray ob- 
servations, therefore we have only a determination of the upper 
limit to their X-ray luminosity. The following analysis is based 
on the sample of 591 low-mass members. 

To compare the luminosity of different classes of objects, 
we analyzed the cumulative X-ray Luminosity Functions (XLF) 
of the Class II and Class III members. In order to take into 
account non-detections, XLFs were d erived using the Kaplan - 
Meier maximum likelihood estimator dKaolan & Meierl Q958). 
Figure [5] shows the comparison of the cumulative X-ray lumi- 
nosity functions of Class II and Class III members. Looking at 
this global description, we note that the range of luminosity is 
3 ■ 10 29 - 3 ■ 10 31 erg s~ l . It is evident that the Class III stars 



8 



M. Caramazza et al.: Coronal properties of NGC 1893 



Table 2. Number of detected X-ray sources and AGN candidates in the 4 Chandra ACIS-I chips. 



CHIP 


Number of X-ray Sources 


Number of AGN Candidates 


10 


158 


28 (0 with flares) 


11 


268 


24 (5 with flares) 


12 


143 


30(1 with flares) 


13 


452 


46 (10 with flares) 



Table 3. Upper limit to the X-ray luminosities for NGC 1893 X-ray undetected members 



ID {Prisinzano et al., 201 1) 
(1) 


a (J2000.0) 
(deg) 
(2) 


<5 (J2000.0) 
(deg) 
(3) 


log L, c 
(ergss ') 
(4) 


408 


80.562999 


33.5 /5532 


30.4/ 


316 


80.566272 


33.449992 


30.31 


454 


80.580917 


33.562509 


30.56 


123 


80.580957 


33.438989 


30.22 


57 


80.590203 


33.447045 


30.29 


213 


80.598221 


33.455573 


30.08 


492 


80.598731 


33.344828 


30.42 


342 


80.602976 


33.532386 


30.24 


497 


80.603157 


33.522741 


30.26 


458 


80.606043 


33.386557 


30.27 


508 


80.613344 


33.350895 


30.35 


134 


80.617675 


33.454839 


30.11 


319 


80.618691 


33.409851 


30.09 


533 


80.620025 


33.347069 


30.61 


439 


80.620694 


33.359367 


30.32 


164 


80.624260 


33.532496 


30.06 


146 


80.624662 


33.559724 


30.16 


267 


80.626424 


33.450228 


29.91 


443 


80.628081 


33.385558 


30.16 


334 


80.637005 


33.457892 


30.09 


418 


80.637332 


33.361037 


30.26 


93 


80.640689 


33.376491 


30.17 


553 


80.641242 


33.370493 


30.2 


452 


80.641437 


33.378881 


30.15 


422 


80.642576 


33.436029 


29.9 


546 


80.644803 


33.350252 


30.67 


243 


80.647512 


33.363512 


30.71 


459 


80.648981 


33.430316 


29.9 


469 


80.650405 


33.388107 


30.06 


522 


80.660049 


33.399277 


30.01 


575 


80.662935 


33.419157 


29.85 


390 


80.663318 


33.413863 


29.8,8 


348 


80.663623 


33.396689 


30.02 


394 


80.671760 


33.378814 


30.28 


337 


80.676571 


33.476855 


29.7 


466 


80.677781 


33.395425 


30.01 


506 


80.679261 


33.379202 


30.12 


136 


80.680467 


33.599258 


30.31 


496 


80.680484 


33.396474 


30.1 


544 


80.681921 


33.495006 


29.7 


296 


80.684036 


33.394367 


30.18 


501 


80.684164 


33.391812 


30.05 


381 


80.684924 


33.393304 


30.02 


524 


80.689466 


33.559017 


30.39 


168 


80.689949 


33.412581 


30.16 


476 


80.690269 


33.350649 


30.23 


472 


80.691751 


33.341261 


30.28 


456 


80.692317 


33.399608 


30.15 


480 


80.693894 


33.367976 


30.13 


54 


80.697590 


33.389405 


31.02 


320 


80.697977 


33.395933 


29.94 


434 


80.698224 


33.40569 


31.37 


488 


80.699277 


33.409286 


29.8,8 


362 


80.701203 


33.354283 


30.21 


165 


80.702890 


33.527785 


30.29 


513 


80.703693 


33.400595 


29.95 


389 


80.704304 


33.416291 


29.82 


427 


80.706539 


33.412613 


29.85 


289 


80.707190 


33.411752 


29.86 


94 


80.707273 


33.428734 


30.37 


413 


80.707867 


33.464616 


29.6,8 


504 


80.708417 


33.437063 


29.71 


401 


80.709488 


33.513873 


29.89 


437 


80.710054 


33.421240 


30.37 


364 


80.710109 


33.462006 


29.71 


447 


80.710368 


33.467710 


29.64 


195 


80.711121 


33.404267 


29.91 


449 


80.711743 


33.4208385 


29.81 


64 


80.712482 


33.567166 


30.09 


59 


80.712499 


33.388352 


29.98 



Note 2. Tablets available in its entirety in the electronic edition. 

Col. (1): Prisinzano et al. ■ Z(J 1 1 ) sequence number, sorted by RA. 

Cols. (2) and (3): Right ascension and declination for epoch (J2000.0). 

Col. (4): Estimated upper limit to the Absorption-corrected X-ray luminosities. 
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Fig. 5. Cumulative X-ray Luminosity Functions of the Class II 
and Class III members. The nondetections were taken onto ac- 
count, deriving the XLFs with the Kaplan-Meier maximum like- 
lihood estimator. 



are brighter in X-rays, the median is higher and the body of the 
Class III object distribution is above that of Class II objects. At 
this age the dynamo mechanism is saturated in most of the PMS 
stars, therefore the X-ray luminosity sca l es as t he bolometric lu- 
minosity of the star (e.g. Preibisch et all l2005h . This means that 
the difference in the whole XLF may be due to intrinsic proper- 
ties or to a different bolometric luminosity (or mass) distribution 
among the two infrared classes. To avoid the mass dependence, 
it is usual to compare the XLFs in several ranges of mass (e.g . 



Flaccomio et 



compare the ALrs in several ranges of mass (e.g. 
all 12003 rjlPreibisch et al J, 12005 ; Prisinzano et al., 



2008), despite the large errors caused by having few objects in 
each bin of mass. We overcome both of these difficulties by 
analysing the Class II and Class III X-ray luminosities as a func- 
tion of the bolometric luminosity of the star. Figure[6]shows the 
quantiles (25%, 50%, 75%) of the Lx distribution in running in- 
tervals of bolometric luminosity covering 80 datapoints, for both 
of Class II (continuous line) and Class III (dashed line) stars. The 
resulting lines have been smoothed over scales of 0.1 in loga- 
rithm of bolometric luminosity. Note that for Class II stars the 
quartiles have been plotted only where the fraction of Lx upper 
limits does not affect the calculation of the quartiles themselves 
(e.g., the first quartile may be calculated only if Lx upper limits 
are less than 25%). The three quartiles give us different infor- 
mation about the distribution of Lx'. the first quartile of Class 
III objects is similar to that correspondent to Class II stars, but 
we are able to calculate this value only for a small fraction of 
Class II member, due to the high fraction of Lx upper limits in 
the faint sample. The median and the third quartile, that describe 
the body and tail of the X-ray distributions respectively, show 
the same feature. They are lower for Class II members but the 
difference decreases with bolometric luminosity and for bright 
stars the quartiles are similar to those of Class III members. 

In order to interpret these results, we also need to compare 
the value of Lx/Lbol for the two infrared classes of objects. We 
show in Fig. [7] the three running quartiles of the distribution of 
Lx/Lbol for Class II and Class III objects. We see that the ra- 
tio decreases for the two classes as a function of Lbol- That 



Fig. 6. X-ray luminosity distribution as a function of the stellar 
bolometric luminosity. The three couples of lines refer to the 
(25%, 50%, 75%) quartiles of the Lx distributions calculated in 
running intervals of Lbol constituted by 80 contiguous points 
and smoothed in interval of bolometric luminosity of 0. 1 erg s" 1 . 
The continuous lines refer to Class II stars, while the dashed ones 
to Class III stars. The quartiles have been plotted just where the 
fraction of upper limit values does not affect the calculation of 
the quartile itself. 



means that the fraction of stars in the saturated dynamo regime 
{Lx/Lbol ~ 3) decreases as a function of Lbol, and the displace- 
ment of Class II Lx/Lbol below that for Class III stars means 
that the fraction of unsaturated X-ray sources is larger for disked 
stars compared to diskless stars. 

We also compared the global X-ray properties of the stars 
in NGC 1893 with those of the Orion Nebula Cluster, in partic- 
ular those obtained from the Chandra Orion Ultradeep Project 
(COUP) dGetman et all 120051) . As we discussed in gf] COUP 
gives the most complete set of studies in X-rays for young stars 
ever achieved for such a rich cluster therefore it has been con- 
sidered a touchstone and its results has been compared with 
previo us and succeeding X-ray studies. iFeigelson & Getrrianl 
(2005) compared the COUP XLF with that of NGC 1333 and 
IC 348 finding that "the shapes of different YSC XLFs ap- 
pear to be remarkable similar to each other, once a richness- 
linked tail of high luminosity O stars is omitted ". Also RCW38, 
IC1396N, NGC 6357 , the Carina Nebula Cluster, M17 and 

i i i i i i i i i ' i 



NGC 2244 (IWolk et all pOOl [German et al 
1 2007 1 : ISanchawala et all 120071: iBroos et 
2008) compared the X-ray properties 



Wang et al 
Wang et al 
of several clusters to 



COUP results, demonstrating that in the Sun neighbourhood the 
X-ray properties of star forming regions are similar, comparing 
clusters of similar age and accounting for the I MF. The only 
releva nt difference in the XLF is for Cep OB3b (Get man et all 
2006), where the XLF has a different shape from that seen in the 
ONC with an excess at Log(Lx) ~ 29.7 or 0.3 M Q . The origin 
of the difference it is not clear but could be ascribable to a de- 
viation in the IMF or some other cause, such as sequential star 
formation generating a non-coeval population. 

Starting from these previous results, we consider the ONC 
as representative of the X-ray properties of nearby star forming 
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Fig. 7. Lx/Lbol distribution as a function of Lbol- The three 
couples of lines refer to the (25%, 50%, 75%) quartiles of the 
Lx/Lbol distributions calculated in running intervals of consti- 
tuted by 80 contiguous points and smoothed in interval of bolo- 
metric luminosity of 0.1 erg s . The continuous line refer to 
Class II stars, while the dashed one to the Class III stars. 



region and will compare COUP results with those of NGC 1893, 
that is ~10 times more distant from the Sun. Figure [8] describes 
the quantiles of the Lx distributions for the Class II and Class 
III stars between 0.35 and 2 M in the case of NGC 1893 (con- 
tinuous line) and ONC (dotted line). The first quantile of the 
NGC 1893 Lx distributions is similar to the ONC, and, since we 
know that the C OUP sample is complete in this range of masses 
dGetman e t al., 2005), it derives that our sample of stars has the 
same completeness of that of the ONC. The median of the two 
Lx distributions are similar, but we note that the third quartile 
is larger for the ONC at the faintest Lbol- This means that we 
lack bright X-ray sources at the lowest masses in NGC 1893 
and this cannot be related to an incompleteness of the sample: 
in that ca se we would loose faint stars an d the effect would be 
enhanced. lAlbacete Colombo et a i] (l2007h . comparing the X-ray 
properties of young stars in Cyg OB2 with those in the ONC, 
noted that dividing the COUP observation in 100 ks segments, 
a significant fraction of the longer and more energetic flares are 
missed. The difference between the two distributions could be 
therefore explained referring to the fact that COUP observation 
is longer than NGC 1893 ones, therefore in the ONC it is pos- 
sible to observe rare very energetic flares that, even not modi- 
fying the median of the distribution, affect its tail. In order to 
take into account the effect of the duration of the observation, in 
the comparison of flares properties between NGC 1893 and the 
ONC (see §3]), we will normalize for the exposure time. 



5. X-ray Variability 

We have analyzed the variability of the lightcurves of the X-ray 
sources, applying the Kolmogorov-Smirnov test to the arrival 
times of photons, and determined that 34% of our 1021 X-ray 
sources (including both members and non-members) are vari- 
able with a probability greater than 95%. If we restrict ourselves 
to only members, the variability fraction is 36 %. When split ac- 




33.e 

Log(L B0L ) [erg s " 



Fig. 8. X-ray luminosity distribution as a function of the bolo- 
metric luminosity. The three pairs of lines refer to the (25%, 
50%, 75%) quartiles of the Lx distributions calculated in running 
intervals of bolometric luminosity constituted by 80 contiguous 
points and smoothed in intervals of bolometric luminosity of 0.1 
erg s _1 . The continuous line refer to the NGC 1893 stars, while 
the dashed one to the ONC stars. 



cording to whether a disk is present of not, 34% of the Class 
III stars are variable vs. 41 % of the Class EE's. While we cannot 
state that these two numbers are statistically different, we have to 
take into account that Class II are intrinsically fainter than Class 
III stars and therefore the comparison has to be done analysing 
sources with the same statistics. Figure [9] shows the number of 
variable stars as a function of the minimum X-ray luminosity 
of the stars in the subsample. We note that, considering just the 
brightest stars, the discrepancy between Class II and Class III 
seems to be larger, suggesting a real difference in the fraction of 
variable stars in the two classes. It is interesting to investigate 
the behavior of Class II stars that still show accretion phenom- 
ena _andJhose_fhatdo notpre sent evidence of strong accretion. In 
the lPrisinzano et all d201 ll) catalog there is a subsample of Class 
II stars for which H„ photometry i s available. Using that tracer 
(see §5.3 in iPrisinzano et all 1201 ll for details on selection), we 
singled out 50 Class II X-ray emitters with ongoing accretion, 
as indicated by H a emission, and 77 that do not show strong H a 
emission. The two samples of Class II X-ray emitters show the 
same fraction of variable light curves. 

A summary of the fraction of variable stars for the different 
subsamples of objects is given in Table|4] 

The Kolmogorov-Smirnov test gives just a global indication 
of the variability of sources, but does not specify the nature of the 
variability. We therefore further analyzed the lightcurves in order 
to study the variability due to stellar flares. To single out flares 
in the lightcurves, we mapp ed the lightcurves using Maximum 
Likel ihood Blocks (MLBs) dWolk et all 120051: ICaramazza et all 
2007). The main characteristic of MLBs is that, being computed 
from the photon arrival times, their temporal length is not based 
on an a priori choice of temporal bin length, but depends on 
the light curve itself; for this reason MLBs are a useful tool 
to quantify different levels of emission, and in particular to de- 
tect short impulsive events, that might be missed by binning the 
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Fig. 9. Number of variable stars as a function of the 
minimum X-ray luminosity of the stars in the sub- 
sample. The vertical bars are the standard deviations 
( VNumber of variable stars/Number of stars in the subsample). 



lightcurves using fixed length bins. Applyi ng the same opera- 
tional definition of flares described in [Caramazza et alj (2007) 
that takes into account both the amplitude and time derivative 
of the count-rate, we singled out flares in the lightcurve as a se- 
quence of blocks with a high count rate and a high rate of vari- 
ation of the photon flux. We estimated the luminosity of each 
flare, scaling the total luminosity of the source during the obser- 
vation by means of a count rate to luminosity conversion factor; 
we then computed the energy of each flare, Efl,, by multiplying 
the flare luminosity for the flare duration, calculated as the total 
temporal length of the blocks associated with the flare. The re- 
sults of this analysis are summarized in Table [4] where we state 
the number of flares per source with energy larger than 2 • 10 35 
erg (this choice will be justified in the following). 

In Fig|T0] we show the number of flares per source as a func- 
tion ofJhe_jrmvimurri X-ray luminosity of the analysed subsam- 
ple. [Caramazza et aQ d2007l) demonstrated that the sensitivity of 
flare detection methods depends on source statistics, therefore it 
is interesting to compare the mean number of flares per source in 
several subsamples of sources, comprising sources with increas- 
ing luminosity. We note that for bright sources where we are able 
to detect all the flares, Class II stars are flaring significantly more 
than Class III stars and therefore that the disk plays an important 
role in generating this kind of high energy events. 

According to the m icroflare s hyp othesis, originally proposed 
for the solar case (see iHudsonl 119911) . the X -ray emission of a 
star can be described as an ensemble of flares with a power law 
energy distribution: 

dN 

=k-EZl with ff>0 (1) 

dE f l r jlr 

where N is the number of flares with energies between Ef/ r and 
Efi r + dEfi r , emitted in a given time interval. If the index of the 
power law (a) is larger than 2 even very high levels of apparently 
quiescent coronal X-ray emission can be obtained from the inte- 
grated effects of many small flares. Figure QT| shows the cumu- 
lative distribution function (continuous line) of the intensity of 
flares for Class II and Class III members. For high counts the dis- 
tribution is well described by a power law, but it progressively 
flattens towards low energies, likely because not all low-count 
flares are individually detectable. 



Fig. 10. Number of flare per source as a function of the minimum 
X-ray luminosity of the stars in the subsample. The vertical bars 
are the standard deviations as in Fig. [9] 

Following the microflares hypothesis, we described the high en- 
ergy part of the differential distribution of flare counts as a power 
law with index a. The cumulative distribution is then described 
by a power law with index a - 1 . 

We have determined the cutoff energy E cu , above which the 
observed distribution is compatible with a power law and the 
relativ e index a — 1 , with the same met hod used bv lStelzer et alj 
(120061) (see also ICrawford et all 1 19701) : 

a - 1 = 1.2 ± 0.2 E cut = 2 ■ 10 35 erg 

In agreement with our assumption on power law shape for the 
distribution, taking E cut larger than the chosen cutoff, the best fit 
value of a remains stable within the uncertainties, while under 
this threshold we cannot neglect the incompleteness effect. 

Figure QT| also shows the cumulative distribution of flare en- 
ergies for the Class II and Class III samples. We noticed that 
Class II star are more variable and more flaring than Class III 
objects, but we now see that the flares of Class II and Class III 
stars show the same distribution. We calculated the slope of the 
two power laws, finding for the two subsamples a- 1 = 1.4+0.4. 

FigureQ~2]shows the comparison between the cumulative dis- 
tribution function of flare energies for NGC 1893 (solid line) and 
the ONC (dotted line). The distributions are normalized to the to- 
tal number of sources of each sample and to the total exposure 
time. Above the completeness cutoff energy the two distributions 
are similar. 

Figure [13] shows the comparison between the mean number 
of flares for the NGC 1893 and the ONC samples. In order to 
compare the same kind of flares, we considered only the flares 
over the completeness cutoff energy. We note that the number 
of flares per source is higher in the ONC but the difference de- 
creases when we exclude low statistic sources. We infer that 
the difference may be due from a likely incompleteness of the 
NGC 1893 sample at low luminosities. 

6. Discussion 

The influence of the environment on star formation process may 
be investigated also through the analysis of X-ray emission. 
NGC 1893, for its position in the Milky Way is a good target to 
investigate possible differences in the behavior of stellar coronae 
and compare stars in clusters at the periphery of the Galaxy with 
those lying in dense spiral arms in the solar neighborhood. 
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Table 4. X-ray variability for different subsamples of X-ray sources 



Sample Number of sources Fraction of variable sources (KS Test 95%) Number of flares per source 

All X-ray sources 1021 0.34 ±0.02 0.16 ±0.02 

Class III X-ray sources 285 0.34 ±0.03 0.19 ±0.02 

Class II X-ray sources 149 0.41 ±0.05 0.27 ±0.04 

Accreting Class II X-ray sources 50 0.40 ±0.09 0.30 ±0.08 

Not accreting Class II X-ray sources 77 0.41 ±0.07 0.27 ±0.06 




Fig. 11. Normalized cumulative distribution function of flare en- 
ergies for Class III (solid line) and Class II stars (dashed line). 
In the range of energies above the cutoff value E cut = 2 ■ 10 35 
erg (vertical segment), the two distributions are compatible with 
a power law with a — 1 = 1.4 + 0.4; for low counts the distri- 
butions flatten, most likely because the detection of low-counts 
flares is incomplete. 



Fig. 12. Cumulative distribution function of flare energies for X- 
ray members in NGC 1893 (solid line) and ONC (dotted line). 
The distributions are normalized to the total number of sources 
of each sample and to the total exposure time. In the range of 
energies above the cutoff value E cut = 2 • 10 35 erg (vertical seg- 
ment), where both the flare energy distributions are complete, 
the distributions of the flares follow similar distributions. 



We analyzed the X-ray properties of NGC 1893, in particu- 
lar comparing the XLF of low mass (< 2M ) Class II and Class 
III young stellar objects. We confirmed the known result that the 
Classical T Tauri Stars (CTTS), i.e. the stars with signatures of 
disks, are globally less active than the Weak-lined T Tauri stars 
(WT TS), i.e. stars that have no signatures of optically thick disks 
(e.g. lPreibisch et aTLl2005h . while being more variable. We fur- 
ther investigated this result, analysing the behavior of Lx as a 
function of the bolometric luminosity. The correlation between 
these variables is compatible with the well known saturated dy- 
namo scenario, typical of pre-main sequence stars. Using the 
method described in S01 we demonstrated that, even considering 
stars with similar Lbol, the median Lx is lower for Class II stars, 
even if the difference decreases for stars with high Lbol- Our sta- 
tistical analysis suggests that the fraction of unsaturated Class II 
members is larger than for the Class III members. This result is 

in agreement with the scenario in which accretors stars are less 

j 1 1 1 

luminous than n o n-accre t ing stars (Stelzer & Neuhauser, 2001; 

Flaccomio et all l2003at ISta ssun et al l 120041 iPreibisch et all 

20051: iFlaccomio et all I2OO6I iTelleschi et all 120071) . However 

the reason why Class II should be underluminous in X-rays is 



not well understood. There are four main ideas that have been 
formulated to explain the different behaviors of disked and disk- 
less stars; two are related to the presence of disk itself and two to 
the presence of accretion from the disk: 1) The lower X-ray lu- 
minosity of Class II stars could be related to the higher extinction 
due to X-ray absorption by circumstellar disks, even if the COUP 
survey results do not support this idea (Preibis ch et all 12005). 
2) A possible explanation is related to the magnetic braking: in- 
deed there are examples of magnetic connection between the star 
and the disk leading to the id ea that disked stars are on aver- 
age s l ower than diskless stars dFavata et all 120051; iRebull et all 
2006; Prisinzano et al., 2008) and have a weaker dynamo ac- 
tion with consequent lower X-ray emission. In this case, the 
dynamo process would be the same for the star with and with- 
out disk and the lower luminosity of CTTS could be attributed 
to the slower rotation of CTTS related to the presence of the 
disk. 3) The presence of accretion may alter the stellar stru cture 
and change the dynamo process itself ( Pr eibisch et all 120051) . in 
this case the dynamo process could be different with different 
satura tion limit for CTTS and WTTS. 4) Finally. iGregorv et alj 
(120071) demonstrated that the corona is not able to heat all the 
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Fig. 13. Mean number of flares per kilosecond above the cutoff 
value E cut — 2 ■ 10 35 erg for subsamples of X-ray sources with Lx 
larger than a given value. Circles refer to the NGC 1893 X-ray 
members, while diamonds to the ONC members. The vertical 
bars are the standard deviations as in Fig. [9] 



material in the accretion columns and this cooler material, not 
visible in X-rays may obscure the line of sight to the star, re- 
ducing the X-ray emission. Recen tly, this scenario was also sup- 
ported by iFlaccomio et all Q2010) who found a correlation be- 
tween the X-ray and optical variability for CTTS, while they did 
not find any correlation for WTTS. With our sample of data, 
we are not able to prefer one of the previous hypothesis, but we 
can add the observational constraint that the difference between 
the two classes of object is more evident for low mass stars and 
we can infer that this can be related to the different evolution 
time of disks around low mass stars and solar mass stars, causing 
there to be a larger fraction of disked stars at the lowest masses 
dCarpenter et all l2006t ISicilia-Aguilar et all 120061 lLada et all 
2006; Megeathetail l2005h . 

We have also compared the X-ray properties of NGC 1893 
with those of the ONC low mass members. At this age the dy- 
namo mechanism is saturated for a large fraction of the stars 
and this leads to a dependence of Lx on the bolometric luminos- 
ity (e.g. lPreibisch et all 12005). Since we cannot a priori assume 
that the mass distribution in the two clusters is similar, we com- 
pared the two samples as a function of the bolometric luminos- 
ity. With this method, we obtained that the median of the two Lx 
distributions vs mass are similar, but the third quartile is larger 
for the ONC at the faintest Lbol- This effect seems to indicate 
an intrinsic lack of bright X-ray sources at the lowest masses 
in NGC 1893. This difference with the ONC cannot be an indi- 
cation of an incompleteness of the NGC 1893 sample because, 
in that case, we would lose faint stars and the effect would be 
in the opposite direction. This difference between the two clus- 
ters in the high luminosity tail of the Lx distributions could be 
ascribed to the longer duration of COUP observation: indeed, 
in the ONC there was the possibility to observe long energetic 
flares that enhance the third quartile of the distribution. 

Analyzing the variability of Class II and Class III stars in 
our sample, we find the known results that Class II are more 
variable than Class III stars and show more flares. The explana- 
tion of this beh a vior c an be related to the effect of accretion. 
IFlaccomio et al ] (1201 Oh observed that the X-ray variability of 
disked stars can be due to the shielding of significant fraction of 
the c oronal plasma by de nse accretion streams of cold material 
(see lGregorv et alll2007l) . It is difficult to compare in detail the 



variability of NGC 1893 sources with the Orion Nebula Cluster 
ones, because of the several biases due to the different distance 
and the different duration of the X-ray observations that lead 
to a different completeness of the samples. Taking into account 
these several biases we find that the X-ray flares in Orion and 
NGC 1893 show a similar energy distribution frequency, leading 
to indistinguishable behavior of the two clusters from the coro- 
nal variability point of view. 

We conclude that despite its peculiar location in the Galaxy, 
NGC 1893 includes a rich population of Class II and Class III 
X-ray sources that have, both from the X-ray luminosity and X- 
ray variability point of view, similar properties to nearby star 
forming regions, such as the ONC. The fact that the X-ray prop- 
erties of clusters in such different environments are so similar 
gives strong support to the Feige lson & Getmanl d2005l) sugges- 
tion that the X-ray properties can be used as standard candles, 
providing a ne w instrument for measurement for distances of 
young clusters dKuhn et aUl2010l) . 

7. Summary and Conclusion 

As part of the large multiwavelenght project The Initial 
Mass Function in the Outer Galaxy: the star forming region 
NGC 1893, we have analyzed the 450 ks Chandra Observations, 
and studied the X-ray properties of the 1021 detected sources. 
The X-ray data have been combined to optical, near and mid- 
infrared data (Prisi nzano et all 1201 lh in order to correlate the 
X-ray properties to the presence of disk and/or accretion. Below, 
we summarize the overall results of our investigations. 

- We derived the X-ray luminosity for our 1021 sources tak- 
ing advantage of spectral fitting for 311 X-ray bright stars 
and from quantile analysis of the other sources, finding a 
range of X-ray luminosity 10 29 5 - 10 315 erg s 1 . We have 
analyzed, in particular, the low mass stars X-ray luminosity 
taking into account the classification of candidate members 
in Class II and Class III m embers based on infrared excesses 
of lPrisinzano et al.ld201 lh . Class III stars appear intrinsically 
more X-ray luminous than Class II stars, even comparing 
stars with the same bolometric luminosity. This effect may 
be due to the presence of the "magnetically connected" disk 
itself or to the ongoing accretion from the disk to the star. 

- We have evaluated the variability of X-ray lightcurve us- 
ing the Kolmogorov-Smirnov test, finding that 34% of the 
sources appear to be variable. We have also searched for 
flares in our lightcurves finding 0.16 flares per source. We 
found that Class II stars are more variable and more flaring 
than Class III stars, while the flare energy properties are the 
same. This property may be related to the accreting cold ma- 
terial that obscures part of the X-ray emitting material, mak- 
ing accreting stars X-ray lightcurves more variable than the 
ones of diskless stars but also suggests that the presence of a 
disk plays a role in generating high energy events. 

- Comparing the X-ray properties of NGC 1893 with those of 
the nearby star forming region ONC, we find that the X-ray 
properties in NGC 1893 are not affected by the environment 
and that a stellar population in the outer Galaxy may have the 
same coronal properties of nearby star forming regions. This 
provides a strong evidence of the universality of the X-ray 
properties, and, as a consequence, a useful tool to determine 
properties, such as the distance or the total population, of 
young clusters. 
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